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Pesticides used improperly can be injurious to human beings, animals, 
and plants. Follow the directions and heed all precautions on labels. Store 
pesticides in orginal containers under lock and key—out of the reach of children 
and animals—and away from food and feed. 


Apply pesticides so that they do not endanger humans, livestock, crops, 
beneficial insects, fish, and wildlife. Do not apply pesticides where there is 
danger of drift when honey bees or other pollinating insects are visiting plants, 
or in ways that may contaminate water or leave illegal residues. 


Avoid prolonged inhalation of pesticide sprays or dusts; wear protective 
clothing and equipment, if specified on the label. 

If your hands become contaminated with a pesticide, do not eat or drink 
until you have washed. In case a pesticide is swallowed or gets in the eyes. 
follow the first aid treatment given on the label, and get prompt medical 
attention. If a pesticide is spilled on your skin or clothing, remove clothing 
immediately and wash skin thoroughly. 


NOTE: Some States have restrictions on the use of 
certain pesticides. Check your State and local regula- 
tions. Also, because registrations of pesticides are under 
constant review by the U.S Environmental Protection 
Agency, consult your local forest pathologist, county 
agriculture agent, or State extension specialist to be sure 
the intended use is still registered. 
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FOREWORD 


Program WIND (Winds In Non-Uniform Domains) was a cooperative 
applied research program conducted by the USDA Forest Service 
(FS) and the U.S. Army, and supported by other public and private 
cooperators. The purpose of Program WIND was to obtain data for 
the evaluation of computer-based predictive mesoscale 
meteorological and particulate dispersion models. The program 
was conducted in northern California over an area approximately 
80 km x 80 km bounded on the south by Sacramento, on the north by 
Redding, on the east by the Sierra Nevada Range foothills, and on 
the west by Interstate 5, and consisted of four phases - Summer 
Phase 1, Winter Phase II, Spring Phase III and Fall Phase IV. 

The field work began in the summer 1985 and ended in the fall 
1987. -The USDA Forest Service participated by conducting 

Siro ate wake and dispersion studies during Phases I, III and IV. 
Miveewas. the first meteorlogical mesoscale study of this scale in 
North America in terms of the type, number, and frequency of 
meteorological data collection. Cooperation among the 
participating agencies and the advancement in meteorological 
instrumentation made this study possible. 


The FS participated actively in Program WIND by providing a 
peegremeciwmector (Jack Barry, Washington Office, Forest Health 


Proreceion start), a public affairs officer (Ann Westling, Tahoe 
National Forest), and assistant program director (Bob Ekblad, 
Missoula Technology Development Center). The FS provided 


administrative management, coordination with the public and news 
Mec vommend contracting and leasing of real estate for the 
numerous study and sampling sites. 


The program approach was to characterize atmospheric conditions 
by obtaining meteorological data during each of valley wind 
patterns associated with the four seasons and during different 
times of the day. Data were used to evaluate and enhance a 
hierarchy of predictive meteorological models. Prior to Program 
WIND the models had suffered for need of a database extensive 
enough to characterize and quantify the physical processes at 
play in the atmosphere. These include models and integrated 
models that predict (1) wind flow over a mountain-valley complex, 
(amine luence of forest and crop canopies on wind flow, (3) 
dispersion of smoke, air pollution, dusts, seeds, pesticides and 
other agricultural materials and (4) effects GieatreraLe wake on 
dispersion of agricultural and forestry sprays and dry materials. 


Spectitcally the PS interest focused on models which predict 
dispersion of smoke from controlled burns and movement Ofea hic 
Doliieants, spray drift, penetration, and deposition of sprays 
into and on forest/range canopies. Data collected during Program 
WIND was also used by the FS to develop and to make FS spray 
dispersion models simpler and easier to use in, Che field. 
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The U.S. Army was primarily interested in a capability to predict 
winds as-they affect obscuration from smoke and dust over an area 
ranging from one km to 6400 km. 


While the FS reported most of its findings from Program WIND at 
symposia, in professional journals, and in FS reports, some data 
escaped reporting. This report and others being issued under the 
Forest Health Technology Enterprise Team (FHTET) and Missoula 
Technology & Development Center (MTDC) report series covers some 
of those unreported studies conducted by the FS in conjunction 
micherProgram WIND. By providing reports within the FHTET and the 
MTDC report series, the data will be more accessable to a broader 
audience through FS and other library searches. 


Key Words: 


forest meteorology, spray models, dispersion models, aerial 
BpOlteacton, forest canopy, spray drift, smoke behavior, forest 
pest management, program WIND. 


ponn W. Barry Harold Thistle 
Director, Forest Health Program Leader 
Technology Enterprise Team - Davis MIDE 


March 1996 March 1996 
















VIS... OAUSS r heswekeuit t2 + reerityy « 
me teve! eu ben. SA ers 161 It “wees. 2028 
i. OARS oF gab 


ATW pattie out *e €tbakt e3<. to-Seoe Bess 
rons brea , Sle vOT Frsrlaedtose 
aie: sset oaded> evan big Date ey @ilet igits ‘ 
hewitt Ba, eehy ; iy 4¢!) ¢ieye ine rey MAGSaT"' 2tee 
woo eolzn¢ 3 nT!) aagee? Sosmeoleved 2's 

i yy) tid J *ein@edw An ene sn boas o geen 


is toe “OTH : Som> |! bee sVvoe we 
“ é ; sl Llinw @9ah. ay 
eri 2 Vrettli tee San £4 


a iitsolt Seno” 
=e, 
, iT os oi. "SS 4 & Ps che ivr “2 + 


Summary 


Evaporation data from Colorado State University have been curve fit to generate 
compatible quadratic relationships for input into USDA Forest Service deposition codes. 
Tables of model coefficients and step-by-step instructions for extrapolating the data are 
included to provide easy program entry for the applicable spray solution. 











Ore 


7 


ca 


~ 
~ 
w 
‘ 
o 
= 
y 
r 
os 
be 
A 


' 
SLETOODD, € : rc Sui ‘vie rviaty matt neering} fi cy ence c 
rshog waters) soetel paved ACIS) Gn wayne We a pint tested SY 
ri inh sof gebelogewes it enolnantue, yon-pegeie Ree meen _ “ + 
DO Yer shies ety ae cond aN avi hy 
| iti 





Data Procedure and Results 


Dennison and Wedding of Colorado State University recently conducted a wind 
tunnel study of evaporation rates of pesticide droplets (Ref. 1). Thirteen pesticide mixtures 
plus water were studied under simulated free fall conditions to extract the evaporation time 
history of the droplets. An extensive test matrix (composed of four droplet diameters, four 
air temperatures and three relative humidity levels) was invoked, with all results stored on 
magnetic tape. An overall error estimate of nine percent was anticipated in the test results 
by Dennison and Wedding due to the test procedure and facility. A portion of the data was 
presented in Ref. 1. A summary of the solutions tested is compiled in Table 1. 


As a consequence of this work, both USDA Forest Service deposition codes, 
AGDISP (Ref. 2) and FSCBG (Ref. 3), have been modified to include the option of 
entering the effects of evaporation on droplet diameter by specifying the three model 
coefficients A , B and C in the equation 


Dan (Q=A+BeeCr 


where 


Diam = droplet diameter (microns) 
t = time (seconds). 


In pursuing the original data on magnetic tape, it was discovered that the desired 
tapes had been overwritten. Thus, only the data plotted in Ref. 1 could be reduced to 
model coefficient form. The digitized data were generated by B. Thompson (private 
communication), and a least squares technique was programmed to determine the A, B 
and C coefficients for best fit of the data. A typical comparison with data is shown in 
Figure 1. The coefficient data are presented in Table 2. 
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11 


ie 
13 
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Table 1: Description of Solutions Tested (Ref. 1) 


Density (gm/cm?) 


0.88062 


1.03551 


0.98770 


0.86090 


1.12600 


1.02004 


1.04280 


1.01224 


1.00882 


1.05404 


1.03008 


1.0 
1.00070 


1.00140 


Mixture 


0.5 gallons of Sevin-4 oil 
0.48 gallons of No. 2 Fuel oil 
0.02 gallons of Automate Red B dye 


0.5 ounces of Polyhedra Douglas Fir Tussock Moth 
Nucleopolyhedrosis Virus 

0.25 gallons of Molasses 

0.7 gallons of water 

NaOH in an amount sufficient to make the pH of the 
mixture between 6.0 and 7.2 


0.8 gallons of Sevin-4 oil 

0.48 gallons of No. 2 Fuel oil 

0.02 gallons of Automate Red B dye 
0. 

0. 

0. 


98 gallons of No. 2 Fuel oil 
02 gallons of Automate Red B dye 


5 gallons of Thuricide 16B 
0.5 gallons of water 


0.5 pounds of Orthene 75S 
1.0 gallons of water 


1.0 pounds of Orthene 75S 
1.0 gallons of water 


2.0 quarts of Dow Esteron (2,4,5,-T) 
10.0 gallons of water 


2.0 quarts of Dow Esteron 99 
10.0 gallons of water 


0.5 pounds of Dipel 

0.5 pounds of Shade 

1.0 gallons of water 

0.016692 pounds of Rhodamine B extra 


0.5 pounds of Dipel 

1.0 gallons of water 

0.016692 pounds of Rhodamine B extra 
water 


5.0 ounces of Nalco-Trol 
100.0 gallons of water 


10.0 ounces of Nalco-Trol 
100.0 gallons of water 


Diameter (microns) 





Time (seconds) 


Figure 1: Comparison of tunnel evaporation data with curve fit for Solution 12; initial 
diameter = 100 microns, air temperature = 20 deg C, relative humidity = 20 
percent. The least squares formula is 

Diam (t) = 100.0 - 6.1375 t - 0.3343 t2 
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(Ref. 1) Index 
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Table 2: Model Coefficients 


Rel Hum 
(percent) 


20.0 
20.0 
20.0 
20.0 


20.0 
20.0 
20.0 
20.0 


20.0 
60.0 


20.0 
20.0 
20.0 
20.0 


20.0 
20.0 
20.0 
20.0 


20.0 
60.0 
90.0 


20.0 
20.0 
20.0 
20.0 


20.0 
20.0 
ZO 
20.0 


20.0 
60.0 


20.0 
20.0 
20.0 
20.0 


A 
(um) 


60.0 
110.0 
240.0 
410.0 


100.0 
100.0 
100.0 
110.0 


110.0 
100.0 


50.0 
130.0 
270.0 
410.0 


110.0 
120.0 
120.0 
130.0 


130.0 
140.0 
130.0 


60.0 
100.0 
240.0 
420.0 


110.0 
100.0 
110.0 
100.0 


100.0 
100.0 


50.0 
110.0 
280.0 
330.0 


B 
(um/sec) 


-0.1147 
-0.0997 
-0.1469 
-0.1131 


-0.0316 
-0.0379 
-0.0471 
-0.0931 


-0.0999 
-0.0863 


-8.7658 
-6.2676 
-4.8895 
-4.0946 


-2.0800 
-3.1428 
-5.2166 
-6.4201 


-6.1130 
-3.7086 
-2.4651 


-0.1033 
-0.0686 
-0.0646 
-0.1598 


-0.0390 
-0.0383 
-0.0507 
-0.0639 


-0.0700 
-0.0835 


-0.0513 
-0.0666 
-0.1664 
-0.1793 


C 
(m/sec) 


0.0005795 
0.0002564 
0.0005760 
-0.0000512 


0.0001074 
0.0000668 
0.0000239 
0.0001818 


0.0002656 
0.0002405 


0.6221 
aly 32 
0.0421 
0.0183 


0.0267 
0.0532 
G71395 
0.1965 


0.1696 
0.0364 
0.0317 


0.000507 1 
0.0001601 
-0.0003182 
-0.0002375 


0.0000679 
0.00007 10 
0.000082 1 
0.0001005 


0.0001761 
0.0001915 


0.0000533 
0.0000990 
0.0006634 
0.0006621 


Figure Solution Air Temp 
(Ref. 1) Index 
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Table 2: Model Coefficients (continued) 


(deg C) 


5.0 
10.0 
15.0 
20.0 


20.0 
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Rel Hum 
(percent) 


20.0 
20.0 
20.0 
20.0 


20.0 
60.0 


20.0 
20.0 
20.0 
20.0 


20.0 
20.0 
20.0 
20.0 


20.0 
60.0 
90.0 


20.0 
20.0 
20.0 
20.0 


20.0 
20.0 
20.0 
20.0 


20.0 
60.0 
90.0 


20.0 
20.0 
20.0 
20.0 


20.0 
20.0 
20.0 
20.0 


A 
(Um) 


100.0 
100.0 

90.0 
110.0 


110.0 
100.0 


50.0 
100.0 
290.0 
380.0 


100.0 
100.0 
100.0 
100.0 


100.0 
100.0 
100.0 


50.0 
150.0 
270.0 
420.0 


130.0 
140.0 
140.0 
150.0 


150.0 
130.0 
120.0 


50.0 
130.0 
240.0 
380.0 


120.0 
130.0 
160.0 
130.0 


B 
(m/sec) 


-0.0146 
-0.0351 
-0.0388 
-0.0636 


-0.0680 
-0.0747 


-4.4632 
-4.4920 
-4.4497 
-4.0137 


-2.6975 
-3.6708 
-4.4740 
-6.2167 


-5.0247 
3.0050 
-1.6839 


10.2330 
-7.4268 
-4.9502 
-3.4388 


-1.8332 
-3.4821 
-4.7848 
=U 959) 


-6.8594 
-4.0980 
-1.4964 


-8.3621 
-6.2391 
-3.1625 
-4.2085 


-2.5886 
-3.1872 
-7.0769 
-5.4953 


C 
(m/sec) 


-0.0000345 
0.000088 1 
0.0000899 
0.000061 1 


0.0001173 
0.0000936 


0.4438 
0.2446 
0.0657 
0.0409 


0.0996 
0.1854 
0.2733 
0.5456 


0.3323 
0.1693 
0.0391 


0.7280 
0.1278 
0.0185 
-0.0029 


-0.01050 
0.00716 
0.03562 
0.03989 


0.09904 
0.03587 
-0.00047 


0.5884 
0.1330 
0.0100 
0.0244 


0.0217 
0.0257 
0.1302 
0.0670 


Figure Solution Air Temp 
(Ref. 1) Index 
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Table 2: Model Coefficients (continued) 


(deg C) 


20.0 
20.0 
20.0 


20.0 
20.0 
20.0 
20.0 
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NWN W NNW Nr 


NNW NRE 


Rel Hum 
(percent) 


20.0 
60.0 
90.0 


20.0 
20.0 
20.0 
20.0 


20.0 
20.0 
20.0 
20.0 


20.0 
60.0 
90.0 


20.0 
20.0 
20.0 
20.0 


20.0 
20.0 
20.0 
20.0 


20.0 
60.0 
90.0 


20.0 
20.0 
20.0 
20.0 


20.0 
20.0 
20.0 
20.0 


20.0 
60.0 
90.0 


A 
(um) 


130.0 
140.0 
110.0 


50.0 
120.0 
310.0 
410.0 


110.0 
110.0 
130.0 
120.0 


120.0 
140.0 
110.0 


50.0 
130.0 
280.0 
400.0 


130.0 
130.0 
130.0 
130.0 


130.0 
120.0 
130.0 


50.0 
130.0 
290.0 
420.0 


110.0 
110.0 
130.0 
130.0 


130.0 
120.0 
110.0 


B 
(m/sec) 


-5.5853 
-5.0716 
-1.9063 


10.0351 
-6.8079 
-5.5876 
-4.5119 


-1.4471 
-3.7666 
-4.8903 
-6.2940 


-4.1637 
-2.7895 
-0.5182 


1225750 
-9.4346 
-4.4417 
-2.1164 


-2.3566 
-2.4665 
-8.9751 
-8.5600 


-7.8133 
-2.6153 
41,9572 


10.5260 
-6.8225 
-6.6731 
-5.0956 


-2.1163 
-4.3289 
-5.8431 
-6.7870 


-6.7514 
-3.9954 
-1.7991 


(m/sec) 


0.0781 
0.0775 
0.0173 


0.6889 
0.1340 
0.0139 
0.0012 


-0.0341 
0.0088 
-0.0023 
0.0855 


-0.0663 
-0.0803 
-0.0245 


LOWS 
0.1851 
-0.0513 
-0.0326 


-0.0502 
-0.0439 
0.1948 
0.1176 


0.0647 
-0.0696 
-0.0570 


09992 
0.2199 
0.0743 
0.0294 


0.0162 
0.0921 
0.1302 
0.2116 


0.2117 
0.0265 
-0.0070 


Figure Solution Air Temp 
(Ref. 1) Index 
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Table 2: Model Coefficients (continued) 


(deg C) 


20.0 
20.0 
20.0 
20.0 


Nr 
SUSY 


NN NN VY Nr NNN WN NN WN 
SSS SUSuA SSSS SSS 


SS SS 


WH NWHHKH NH 
eee CT OOM lol Boole lo Mo lolo lo Boole ol oloro 


SSS SS SS SSS 


NNW Ne 


Rel Hum 
(percent) 


20.0 
20.0 
20.0 
20.0 


20.0 
20.0 
20.0 
20.0 


20.0 
60.0 
90.0 


20.0 
20.0 
20.0 
20.0 


20.0 
20.0 
20.0 
20.0 


20.0 
60.0 
90.0 


20.0 
20.0 


20.0 
20.0 


20.0 
60.0 
90.0 


20.0 
20.0 


20.0 
20.0 


20.0 
60.0 
90.0 


A 
(um) 


50.0 
140.0 
280.0 
390.0 


120.0 
130.0 
130.0 
140.0 


140.0 
140.0 
110.0 


50.0 
100.0 
240.0 
380.0 


110.0 
100.0 
100.0 
100.0 


100.0 
110.0 
100.0 


110.0 
400.0 


120.0 
120.0 


120.0 
110.0 
100.0 


130.0 
400.0 


130.0 
130.0 


130.0 
120.0 
100.0 


B 
(\1m/sec) 


-10.8030 
-7.4536 
-4.3828 
-3.2595 


-2.5664 
-4.0313 
-3.8285 
-6.0584 


-6.1061 
-4.9250 
-2.5407 


-10.8690 
-6.1375 
-3.5386 
-2.5431 


-1.6399 
-3.0005 
-3.5431 
-5.6546 


-5.9072 
-2.3116 
-1.3256 


-5.6661 
-2.1162 


-0.9795 
-3.9073 


-4.0762_ 


-1.1747 
-0.4519 


-4.8199 
-3.2936 


-1.8616 
-3.4725 


-4.0621 
-2.3141 
-0.5218 


(um/sec2) 


0.9753 
0.1724 
0.0237 
0.0038 


0.0222 
0.0426 
0.0122 
0.0596 


0.0675 
0.0436 
0.0275 


-0.9249 
-0.3343 
-0.0259 
-0.0300 


-0.0751 
-0.1750 
-0.2503 
-0.4011 


-0.3439 
-0.1778 
-0.0180 


-0.3342 
-0.0219 


-0.4107 
-0.3326 


-0.3138 
-0.2642 
-0.0593 


-0.2681 
-0.0255 


-0.3361 
-0.4495 


-0.3499 
-0.2712 
-0.0525 


Data Correlation and Extrapolation 


When attempting to make use of the evaporation data included here, the following 
information must be determined: 


TEP= 


The evaporating solution must be consistent with one of the fourteen tested 
solutions, or approximated by one of the fourteen solutions. If the desired solution does 
not fit any of the fourteen possibilities, this document cannot be used to quantify its 
evaporation. 


STEP 2 


If the initial droplet diameter, air temperature and relative humidity are within nine 
percent of any of the tested conditions for the desired solution, the model coefficients can 
be used directly as tabulated in Table 2. If not, a computation must be undertaken to 
determine a set of correlated model coefficients. 


The original experimental data set of Dennison and Wedding contained all 
intersection points in the three-dimensional test matrix cube of 50, 100, 250 and 400 
microns (nominal) droplet diameters; 5, 10, 15 and 20 deg C air temperature; and 20, 60 
and 90 percent relative humidity, for each solution. Even with some dropout points (where 
data were not obtained), they accumulated 496 data runs. This data would have been 
sufficient to interpolate for any tested solution at any initial droplet diameter, air temperature 
and relative humidity. 


Unfortunately, only 143 data sets remain (the data plotted in Ref. 1). These data 
are restricted to 20 deg C and 20 percent relative humidity for all droplet diameters; 100 
microns (nominal) and 20 deg C for all percent relative humidity; and 100 microns 
(nominal) and 20 percent relative humidity for all air temperatures, for each solution. In 
Statistical terms the recovered data display no interactive effects; they lie upon three 
intersecting edges of the test matrix cube. 


To overcome this problem and populate the test matrix, approximations must be 
made that will generalize the B and C coefficients for any diameter, air temperature and 
relative humidity, for the fourteen solutions tested. The procedure is as follows. 


STEPS 
The A model coefficient is the desired initial droplet diameter. 


The B model coefficient is obtained by curve fitting the B coefficients determined 
from test data, and assuming separation of variables. Since the air temperature and relative 
humidity effects are known only for 100 microns (nominal), it has been assumed that their 
effect is similar across all droplet diameters. Thus, if a doubling of the air temperature 
halves the droplet diameter for 100 microns, it will do the same thing for any other initial 
droplet diameter. Evaporation is a complex theoretical problem even for Solution 12 
(water), and is more so with additives. Since evaporation reduces as the air temperature 
goes to zero and the relative humidity approaches -100 percent, a separation of variables 
effect seems the easiest (and most obvious) to impose on the data. Thus the B correlated 
model coefficient is approximated by the expression 


b d 


B=aD°-T°H 


where 


D = initial droplet diameter (microns) 
T = air temperature (deg C) 
H = relative humidity (0 - 100). 


The B model coefficients in Table 2 are plotted logarithmically to recover the least 
Squares power law for each variation. Typical correlation plots are shown in Figure 2. 
Consistent tabulated results are given in Table 3. 


The C model coefficient is more difficult to quantify because it is more difficult to 
extract accurately from the data. Separation of variables will again be assumed. In this 
case the air temperature and relative humidity effects will be taken as linear 


C=eD' 


T \{100-H 
20 }\~_ 80 


The C model coefficients in Table 2 are then plotted logarithmically to recover the 
least squares power law for initial droplet diameter. A typical correlation plot is shown in 
Figure 3. Consistent tabulated results are also given in Table 3. 

EXAMPLE 

For Solution 12 (water) with an initial droplet diameter of 200 microns, air 


temperature of 12 deg C and relative humidity of 50 percent, substitution of the appropriate 
correlation model coefficients from Table 3 gives 


Diam (t) = 200.0 - 1.052 t - 0.028 t2 


as the equation to represent evaporation effects for these specified conditions. 
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Figure 2a: Correlation of B model coefficients for Solution 12 for initial droplet 
diameter. 
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Figure 2b: Correlation of B model coefficients for Solution 12 for air temperature. 
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Figure 2c: Correlation of B model coefficients for Solution 12 for relative humidity. 
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Table 3: Correlated Model Coefficients 


Solution 

Index a b 
1 --0.00719 0.1015 
Dee 660m 620.3577 
3 0.00370 —0.2390 
4 -0.00006 0.7506 
5 -11.342 —--0.0837 
6 -62.021  —-0.5001 
7 32.849 —--0.4233 
878.219) -0.3092 
9 -379.67 0.8332 

10 -23.888 ~—--0.2871 

11 -71.441 —--0.5725 

12 -172.58 0.6829 

(emeit4438 -016243 

14 10.547 —--0.1908 

b 


Diam (t)=D+aD°T°H 


where 


Cc 
0.8458 
0.8670 
0.4525 
0.9991 
0.4758 
0.9854 
0.6678 
0.9603 
1.0691 
0.8613 
0.7124 
0.8716 
2.1690 
ANS Ys) 


t+ 


d 
-0.1031 
-0.5701 

0.1941 

0.1777 
-0.6312 
-0.9070 
-0.5740 
-1.3011 
-1.0303 
-0.7910 
-0.5780 
=0.932) 
-1.4883 
-1.1394 


d.. eD'T(100-H) 2 
1600 


e 
0.0139 
637.54 
0.0002 
0.0 
DUS sey 
20530.0 
678.54 
45644.0 
314.7 
508.98 
23465.0 
-2604.2 
-10550.0 
-28629.0 


Diam = droplet diameter (microns) time history variation 
D = initial droplet diameter (microns) 
T = air temperature (deg C) 


H = relative humidity (0 - 
t = time (seconds) 


100) 
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-0.8166 
-1.7104 

0.0076 

1.5380 
-1.3064 
PRON ST| 
-1.8452 
-2.7627 
-1.5743 
-1.5928 
-2.5379 
aod 29 
-2.1824 
-2.3251 
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Figure 3: Correlation of C model coefficients for Solution 12 for initial droplet 


diameter. 
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